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ABSTRACT 

We examine how coupled dark matter and dark energy modify the development of Zel'dovich pan- 
cakes. We study how the various effects of these theories, such as a fifth force in the dark sector 
and a modified particle Hubble drag, produce variations in the redshifts of caustic formation and the 
present-day density profiles of pancakes. We compare our results in direct simulation to a perturba- 
tion theory approach for the dark energy scalar field. We determine the range of initial scalar field 
amplitudes for which perturbation theory is accurate in describing the development of the pancakes. 
Notably, we find that perturbative methods which neglect kinetic terms in the scalar field equation of 
motion are not valid for arbitrarily small perturbations. We also examine whether models that have 
been tuned to match the constraints of current observations can produce new observable effects in the 
nonlinear structure of pancakes. Our results suggest that a fully realistic three-dimensional simulation 
will produce significant new observable features, such as modifications to the mass function and halo 
radial density profile shapes, that can be used to distinguish these models from standard concordance 
cosmology and from each other. 

Subject headings: cosmology:theory, dark matter, dark energy, structure formation, methods: N-body 
simulations 
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1. INTRODUCTION 

Dark energy is perhaps the most profound and essen- 
tial mystery in modern cosmology. While the ACDM 
cosmological model has proven very successful in explain- 
ing and predicting many features of our universe, such as 
the fluctuations in the co smic microwave background (eg. 
Ide Bernardis et al] 120001) . the l arge-scale matter distri- 
bution (eg. iPercival et al.l l2001h , and distance measure- 
ments to type l a supernovae (jPerlmutter et al.l Il999t 
iRiess et aLlll998D . the physics of the dark sector (dark 
matter and dark energy), which comprises roughly 
ninety-six percent of the energy density of the universe, 
is largely unknown. 

Currently, there are too few observational constraints 
to determine the precise nature of the dark energy and 
any possible int eractions it might have with dark and 
baryonic matter (|Bean et al.ii,20 05l. However, we can use 
simulations of nonlinear structure formation to explore 
the consequences of plausible dark energy models, includ- 
ing those that propose couplings between dar k matter 
(DM) and dark energy (DE) fsee lAlcanij|2006l and ref- 
erences therein). Models that have pa rticle physics mo - 
tivations often predict such couphngs () Amendolal l2000f ) . 
Such theories are intriguing because they might pro- 
vide a resolution to current cosmologi cal problems, such 
as th e so-c alled coincidence problem ( Zimda hl fc PavonI 
I200H lAmend ola fc Tocchini-Valcntini 20 01), and the 
obser ved emptiness of the voids (jFarrar fc Peebles! 
l2004l), the latter of which was con firmed numerically 
bv lNusser. Gubser. fc Peebles! (|2005( ). 
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An important goal of coupled DM-DE simulations is 
to identify observational methods to test these theo- 
ries. The observables studied to date include luminosity 
distances (lAmendola et al.!!2007f). t he growth of rn atter 
perturbations ( Olivares et al.l2006l; iK oivisto |2005! ). the 
abundance of clust ers dM anera fc Mota 2006), and the 
Sandage-Loeb test (jCora saniti et al. 2007). Not only can 
we use such simulations to search for additional observ- 
able features, but we can also use them to find ways to 
distinguish models and determine the validity of pertur- 
bation methods. 

Much of the work to date has follow ed the framework 
established bv lFarrar fc Peeb les (2004). This model uses 
a dynamical scalar field to provide the dark energy, and 
it allows that field to couple to dark matter via a Yukawa 
interaction. Although there are som e issues with models 
of this type (jDoran fc Jackel!!2002f ). they provide useful 
foils for studying the possible effects that other mod- 
els, which are complicated but more robust, might pre- 
dic t. These models would include the model discussed 
by !Huey fc Wandelt! (|2006! ) and the two-family model 
considered in Farrar and Peebles. 

Some of the previous work has examined structure 
formation with modifications due to a fi fth force, al- 
beit at low spatial and force resolution (|Nusser et al.! 
12005!: iMaccio et all I2004D . However, we believe that 
an incremental approach that analyzes the various ef- 
fects of DM-DE interactions on simpler initial condi- 
tions is crucial before tackling more realistic initial con- 
ditions. This approach allows the effects observed in 
more realistic simulations to be understood and gener- 
alized. Therefore, in this paper, we will study the ef- 
fects of coupling between dark matter and dark energy 
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on the development of Zel'dovich pancakes. Zel'dovich 
pancakes are the well-known solutions to the problem of 
the gravitational collapse of one-dimensional, sin usoidal, 
plane - wave, adiabatic density perturbations (Zer dovichI 
Since they are well-studied in a variety of contexts 
not i nvolving DM-DE interactions (for example. .Morei 
IQSaiGnedinlllQQaiYuan et allflOQlt iValinia et aLl fTooa 



Anninos fc Normanlll994l; lAnninos et aLlllOOsT and oth- 
ers), we can more easily understand the effects of addi- 
tional physics on structure formation, laying the ground- 
work for a more complete three-dimensional study. 

In Section [2] we discuss the relevant equations, the ef- 
fects we will study, and our numerical techniques. Sec- 
tion [3] discusses the role that the exotic physics intro- 
duced by the Farrar and Peebles model, such as a DM 
particle fifth force and time-dependent DM particle mass, 
plays in structure formation. We are careful to separate 
the various effects predicted by the model. Even if this 
particular example does not withstand closer scrutiny, 
unrelated theories may predict similar effects. 

In addition, so far many researchers have exploited per- 
turbation theory to treat fluctuations in the DE scalar 
field. However, we do not know a priori the validity of 
this approach, especially as we begin to explore the non- 
linear consequences of this theory. In Section |4] we will 
compare the structure formation results from perturba- 
tion theory alone to the results from doing a complete 
nonlinear analysis. We will also discuss the regime where 
perturbation theory is most valid in accurately predict- 
ing structure. 

Many of the theories of this type have adjustable pa- 
rameters, and these parameters must be adjusted to 
match ACDM predictions, at the risk of being indistin- 
guishable from it. There are also usually several unique 
combinations of parameters that provide similar, if not 
identical, results. In Section[5]we attempt to find ways to 
distinguish models that remain indistinct in perturbation 
theory. Also, we will determine if nonlinear effects can 
distinguish models from standard cosmology even when 
effects based on perturbative methods cannot. 



2. THE MODEL 



2.1. 



Coupled Dark Matter and Dark Energy 

Following Farrar and Peebles, we will consider a dark 
energy (DE) scalar field (j) with action 



S 



DE 



]^<t>,,r - v{4>) 



(1) 



and a single classical nonrelativistic dark matter (DM) 
particle family with action for the ith particle 



S 



DMi 



y\(t){x,)- (t)^\^Jgf,^dx'^dx'^, (2) 

where y is the diniensionless Yukawa interaction strength 
and 0, is a constant providing an intrinsic mass to the 
DM particle. We will set 0* = 0, so that the DM particle 
mass is due entirely to the field value. In both of the 
above actions and throughout, we have set fi. = c = 1, 
and we will ignore baryons. 

Assuming a spatially flat Fricdmann universe, we can 
use equation Q to obtain a comoving DM particle equa- 
tion of motion in one dimension: 




dx 



4> dx 



(3) 



Here, $ is the normal comoving gravitational potential, a 
is the scale factor, and v is the comoving particle velocity, 
and X is the comoving position. Throughout, dots refer 
to derivatives with respect to proper time t. The field 
action in equation ([1]) gives the evolution of </>: 
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where p is the DM particle comoving density. 
The comoving potential satisfies the Poisson equation: 

V2$ = — p-p. 5 

Here and throughout, an over line indicates a spatial av- 
erage. 

The Friedmann equation, neglecting radiation, curva- 
ture, and baryonic terms, is now 

^ — s 2 



a ' ' 



(6) 



A subscript of denotes the present-day value. The first 
term on the right-hand side reflects the contribution of 
the DM with its time-dependent mass. The terms in 
brackets are, respectively, the kinetic and potential ener- 
gies of the DE scalar field. 

We do not know the initial conditions of the field, but 
we do know that today the field behaves as a cosmological 
constant, so the potential term dominates and has a value 



V{4>o) = flAPcrit- 



(7) 



Also, at early enough times, Farrar and Peebles found 
that equation ([4]) reduces to 
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dt 
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which we use to set the initial condition for 0. 

We can identify four unique ways in which the extra 
interactions modify structure formation. First, the DM 
particle mass directly depends on the field value, so we 
may write the ratio of the modified mass to its present- 
day value as 

moM 4> /„x 
T] = = — (9) 

Secondly, the interactions modify the Hubble drag found 
in the DM equation of motion, so that its ratio to the 
drag in standard cosmology is 
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Next, we notice a modified particle acceleration due to a 
fifth force in equation ([3]). We will define 
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dx dx 

dx 



(11) 



Finally, the dynamic scalar field itself indirectly affects 
structure formation via a time-varying Oa in the Fried- 
mann equation. We define 



(12) 
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If the fluctuations in the field are small enough, we 
may do perturbation theory. Farrar and Peebles found 
that in this regime, we may replace </>(x) with a single 
spatial average, (pb and a sufficiently small perturbation 
field (j)i{x). We may then write equation ([4]) as 

+ 3-(/)b + — + -— — -^a = 0. 13) 

a a(pb STTUcpbfi 

Also, the fifth force ratio appears instead as 

+ (14) 

We will contrast our results with standard concordance 
cosmology, which we achieve by setting to a single value 
satisfying equation ([7]) and by preventing from evolving 
dynamically. 

We have the freedom to choose an appropriate poten- 
tial V{4>). Although there are m any potentials in the lit - 
erature, suc h as the exponential (|Ratra fc Peebleslll988h . 
power- law (Peebles & Ratra"1988'), and power-law and 
sine (Dodclson ct al. 2000), we will adopt the power-law 
potential found in Farrar and Peebles: 

y(0) = K/r, (15) 

where we are also free to choose the constants K and 
a. Potentials like this lead to reasonable behavior, 
i.e. potential-dominated solutions at the present epoch. 
Again, we chose this potential merely as an example. 

2.2. Numerical Techniques 

For our simulations we developed a one-dimensional 
iV-body particle-mesh code. We used cloud-in-cell 
mappi ng for interpolating between the mesh and par- 
ticles ()Hocknev fc EastwoodI I1988D . and a second or- 
der leapfrog integration scheme for particle advance- 
ment. In one dimension and with finite differenc- 
ing, we can solve the Poisson equation exactly using 
a Thomas algorithm (|Conte fc dcBoor 1972), modified 
for pe riodic boundaries via the Sherman-Morrison for- 
mula ([Press et al.lll992f l. We discuss the details of our 
scheme for solving the scalar field and scale factor in the 
appendix. 

For all calculations, we used ilm = 0-26, Ha — 0.74, 
and Ho = 100 /i = 71 km s~^ Mpc"\ All runs took 
place in a one-dimensional box of length 10 Mpc per 
side, with 65, 536 particles to represent the dark matter 
and 8, 192 zones for the Poisson solver and the scalar 
field. Since we used finite-differencing when forming the 
gradient to construct the fifth force, we required a large 
particle-zone ratio to dampen noise in the density field, 
which couples to the scalar field. 

All simulations used the same initial conditions. We 
distributed particles evenly throughout the grid (de- 
scribed by position qi = iAx) and perturbed them using 
the Zcl'dovich approximation: 

2 

x = q-\ — aAsm(kq) (16) 
5 

2 

v = -aAs'm(kq), (17) 
5 

where A: = 27r/A is the comoving wavenumber of the per- 
turbation. The amplitude A can be written in terms of 
the redshift Zc of the formation of the first caustic: 

A^ + zc)/{2k). (18) 
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-2 


0.03 


1.89 


1.72 


B 


-2 


0.0057 


4.02 


3.94 


C 


6 


2.0 


1.80 


1.68 


D 


6 


280.0 


3.89 


3.83 



TABLE 1 

Simulation parameter choices. 



We chose the initial perturbation amplitude such that 
the first caustic would form at = 5.0 for an fl^ = 1.0 
universe, and we chose the comoving wavelength of per- 
turbations to be A = 10 Mpc. These choices are 
arbitrary, but commonly used in the literature. All com- 
putations start at a redshift of z = 50 with DM particle 
masses determined by equation ([5]). 

We do not know the initial conditions of (j) in advance, 
so we must guess initial conditions and iterate until we 
meet the condition d = Hq at the present epoch. Based 
on the comments made by Farrar and Peebles, we chose 
four combinations of the potential parameters K and a 
that yield reasonable behavior. Table [T] lists the param- 
eters, the guessed initial field value at z = 50, and the 
field value today as calculated from equation (|4]). We 
performed these calculations in perturbation theory. We 
selected these parameters for behaviors that were con- 
sistent with current observations, but provided unique 
evolutions of the scalar field for comparison. 

3. ANALYSIS OF EFFECTS IN PERTURBATION THEORY 

We now illustrate the effects of these additional inter- 
actions on structure formation. To highlight each indi- 
vidual consequence of the DM-DE interaction, we sepa- 
rately analyze the resulting structure when including the 
effects of the variable mass, the modified Hubble drag, 
the fifth force, and the dynamical field. Finally, we will 
examine the final structure when these effects are com- 
bined. In all that follows, we highlight the behavior when 
using the potential parameters K = 0.03 and a = —2 
(parameter set A) as an example. 

Figure [1] shows the evolution of 77, 7, (3pert, and S as 
functions of scale factor for parameter set A. On aver- 
age, the fifth force is about three percent of the gravita- 
tional force. We expect the fifth force to play a significant 
role throughout cosmic history. The Hubble drag is re- 
duced by about one percent. The particle mass starts at 
ten percent greater than its present-day value, but by a 
scale factor of 0.5 it is only two percent larger. We ex- 
pect the modified mass to contribute to structure mostly 
in the early universe, but by the present day its effects 
will not be significant. In the early universe, the dynamic 
field gives a large Q\, and hence a more rapid expansion, 
which should delay caustic formation relative to standard 
cosmology. Notice that our calculations do not give pre- 
cisely the required values of f^A, but are less than two 
percent off. We fixed our field at late times such that the 
energy density of the field was purely in the potential, 
when in reality the kinetic term, while very small, is not 
entirely negligible. Fortunately, this small error does not 
in any way significantly alter structure formation. 

We show the density profile at z = in Figure [H mod- 
ified by the individual additional effects of the DM-DE 
interaction. We do not show the result from an evolv- 
ing mass and dynamical field, because these profiles are 
hard to distinguish from standard cosmology. When iso- 
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Fig. 1. — Comparison of additional effects due to the DM-DE 
interactions. Shown are relative values defined above as rj (solid 
line), 7 (line with hatches), Ppert (line with stars), and 5 (line with 
boxes). The symbols are only to aid in distinguishing the lines. 
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Fig. 2. — Dark matter density profile at z = 0. The solid line in 
all plots is the standard cosmology. The dotted line indicates the 
effects of the additional interactions: the top panel shows only the 
result with the drag term included, and the bottom panel shows 
only the fifth force included. Here, x is the distance from the 
midplane of each pancake. The numerical shot noise in the density 
near x = 1 is negligible and unimportant. 

lating the individual effects of the drag, fifth force, and 
dynamic field, we forced the scalar field to evolve as if the 
scale factor obeyed the modified Friedmann equation in 
equation ([6]); however, the actual scale factor evolved ac- 
cording to concordance cosmology in these comparisons. 
We did this because without the dynamic mass in the 
Friedmann equation, the scalar field would not reach its 
required value today. 

We can easily understand the modifications due to the 
fifth force: an additional attractive force causes the caus- 
tics to form earlier, and hence each peak in the density 
profile at z = is farther away from the center. The 
effect on each peak is similar: a fifth force only in the 
early universe will affect mostly the first two caustics, 
but since the force ratio remains roughly constant even 
in late times, the latest peak also shifts. The overall size 
of the pancake is then larger. 

Due to the modified drag term alone, the first two caus- 
tics form earlier, but the latest caustic forms at the same 
time as in standard cosmology. This requires more ex- 
planation. The Hubble drag term is most important in 
the early universe. Due to the reduced drag, particles 
move faster, and reach a greater turnaround radius. So 
the first two caustics form earlier for the modified cos- 
mology. However, at later times, the drag is no longer 



significant, and the third caustic forms at nearly the same 
redshift as in standard cosmology. Also, it appears that 
modifications to the drag are very important: even small 
deviations in the redshift of caustic formation lead to 
significant difference in the final peak locations. 

The increased mass creates tension between two op- 
posing influences on structure formation. Since the par- 
ticles already in halos are more massive, they create a 
deeper potential well and encourage other particles to 
fall in faster, thereby causing caustics to form earlier. On 
the other hand, a larger in the Friedmann equation 
accelerates the expansion in the early universe, which 
increases the Hubble drag and dampens structure for- 
mation. We found that the expansion effect dominates, 
but not significantly, and caustics form only slightly later 
than in standard cosmology. Also, since the DM parti- 
cle is more massive in the past, the regions near early- 
forming caustics have a higher density, but by 2; = the 
mass is the same as in standard cosmology, so the overall 
amplitude of the density profile does not change. 

We must contrast this evolving-mass universe with 
a ACDM cosmology having simply a larger ^„i- In 
such a universe, the competing effects of increased mass 
roughly cancel each other out, and caustics form at simi- 
lar redshifts to those in standard ACDM. In an evolving- 
mass universe, however, as the mass decreases, the halo 
potential wells get shallower, but the universe has al- 
ready expanded more, so outlying particles take longer to 
reach the halo, and caustics correspondingly form later. 
At lower redshifts, the particle masses are nearly their 
present-day values, so this effect is not as significant, and 
the final peak locations are slightly closer to the pancake 
midplane. 

Finally, the dynamic field has the expected result: 
more rapid early-universe expansion pulls particles away 
from each other, and caustics take longer to form. As ex- 
pected, this effect is nearly negligible in the late universe, 
as the field approaches its present-day value. Unlike 
the differences between an evolving-mass universe and a 
constant-mass universe, the structure modifications from 
a dynamic field are, in this case, nearly indistinguishable 
from a universe with a roughly ten percent larger, but 
still constant, J^a- 

We can also examine the effects on particle velocities 
due to these additional interactions. Figure [3] is a DM 
particle phase plot at z = with modifications due to 
the reduced Hubble drag and the fifth force, compared 
to standard cosmology. The dynamic field and variable 
mass terms have no effect, but as expected, the reduced 
drag and fifth force increase particle velocities. 

When we combine results, in Figure [H we see that 
the effects of the modified mass and the dynamic field 
tend to cancel out the fifth force. So, the final density 
profile most closely resembles the results from the drag 
term alone. The innermost caustic is the same distance 
from the midplane of the pancake, but the outermost and 
middle peaks are roughly 15% farther away. The slope 
of the density profile remains undisturbed and the total 
number of peaks formed by the present epoch are the 
same. 

4. THE VALIDITY OF PERTURBATION THEORY 

So far, we have only considered perturbation theory, 
relying on a single value to describe the DE scalar field 
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Fig. 3. — Dark matter particle phase plots at 2 = 0. The solid 
line in all plots is from concordance cosmology, and the dotted 
lines with open squares are the results from including the effects 
indicated in each plot. The open squares only represent select data 
points, to aid in distinguishing the lines. We have omitted the plots 
with dynamic field and dynamic mass effects, since these have no 
discernible influence. 
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Fig. 4. — Dark matter density profile at z = 0. The solid line 
is from standard cosmology. The dashed line is the density from 
using perturbation theory, and the dotted line is from using full 
theory. This is including all effects. Note that perturbation theory 
is nearly indistinct from the full theory. The numerical shot noise 
in the density near a; = 1 is negligible and unimportant. 

at a specific time. When we relax this constraint, each 
particle feels its own force, has its own drag, and has an 
independent mass. When solving the Friedmann equa- 
tion, we must take a spatial average of the field, rather 
than using a single background value. Also, the expres- 
sion for the fifth force changes to be directly proportional 
to gradients in the scalar field. 

There is some ambiguity when considering the initial- 
ization of the scalar field. While the average value must 
match the perturbation theory background value, we are 
left with few clues to the initial wavelength and ampli- 
tude of the perturbations. So, using adiabatic reasoning, 
we initialized (f) to have the same spectrum as the den- 
sity field, which in our case is a single perturbation mode 
with wavelength A = 10h~^ Mpc at 2: = 0. However, we 
freely chose the amplitude, and we found that to make 
an attractive fifth force, which is required for a Yukawa- 
type interaction, we needed (p to have a phase opposite 
to that of the density field. 

We must take care in the initialization of both the 
background value and amplitude of perturbations of the 
field when considering the full theory. First, the addi- 
tion of the gradient term in the equation of motion for 
the field will affect the field's background evolution, and 
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Fig. 5. — Average ratio of fifth force to gravitational force — 
1) for perturbation theory (solid line), compared to the full result 
with varying initial amplitudes (dotted lines with various marks, 
labeled in figure). The marks on the lines only represent select data 
points, to aid in distinguishing the lines. Shown are the average 
strengths up to the given scale factor. 

may force us to modify the initial value, so that still 
satisfies the present-day condition in equation For- 
tunately, the gradients are not large, and this term does 
not modify the evolution of greatly. Hence, we con- 
structed the field so its average value matched the back- 
ground perturbation value. It appears that perturbation 
theory is perfectly appropriate for probing the linear re- 
sults of this model, such as estimates of the Hubble time 
and the location of peaks in the CMB power spectrum, 
which Farrar and Peebles discuss. 

Secondly, by varying the initial amplitude, we can 
chang e the strength of the fi ft h forc e at high red- 
shift. iKesden fc Kamionkowskil (|2006ai rbl) have shown 
that a constraint on any dark matter fifth force to within 
only a few percent of the gravitational force is obser- 
vationally possible. The fifth force under perturbation 
theory agrees with this constraint: Ppert ~ 1-03 through- 
out cosmic history. Figure [5] shows the space-averaged 
and time-averaged value of Ppert — 1 and /3 — 1, which 
we denote as (/3) — 1, versus redshift for varying initial 
amplitudes. This is the average fraction of gravitational 
acceleration experienced by the DM particles due to the 
fifth force up to the given scale factor. The field fluctu- 
ates rapidly in the early universe, and by time-averaging 
we can better examine the general force trends. The 
fluctuations in the fleld decay rapidly and asymptotically 
approach a constant value by redshift z = 45. 

No matter the initial amplitude, the average perturba- 
tion theory fifth force ratio agrees to within ten percent 
with the full result by the present epoch. However, the 
strength of the fifth force matters much more at high red- 
shifts: different forces here can greatly affect structure 
formation. We found that an initial amplitude contrast 
of S^{z = 50) = 10^* provides roughly the same order of 
magnitude additional acceleration as perturbation the- 
ory predicts at early times. 

For amplitudes smaller than 10~*, the results for the 
fifth force do not change. Since the scalar field couples 
directly to the dark matter density, the field will always 
contain perturbations. If we initialize the amplitudes be- 
low the threshold of 10~^, the density coupling will pro- 
duce fiuctuations in the scalar field spectrum, creating a 
fifth force. However, larger amplitudes provide too much 
of an acceleration at early times. For example, a contrast 
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of d^{z — 50) = 10^® provides an average fifth force that 
is almost half as strong as gravity at early times. We 
found that if the initial contrast is larger than can be 
handled adequately by perturbation theory, then we get 
radically different structure formation, such as a very 
early first caustic and the development of four caustics 
by the present epoch. Note that these results are largely 
redshift-independent: the DM coupling will always end 
up dominating the spectrum of scalar field fluctuations 
after a short period of time. The large fifth force for am- 
plitudes above 10~^ is due to the large artificial initial 
scalar field amplitude. 

However, if the amplitude of perturbations is small 
in the early universe, then perturbation theory undcr- 
emphasizes the fifth force at high redshift. The fact that 
perturbation theory cannot describe large fluctuations 
makes sense. However, we would expect that perturba- 
tion theory should accurately capture the behavior of all 
initial amplitudes below a certain value. 

With the above considerations, we examined the case 
of an initial contrast of S^{z = 50) = 10^®, since this 
provides a fifth force large enough to affect structure, 
but which is not ruled out by observations. We compare 
perturbation theory to the results from a full dynami- 
cal field. We found that for this amplitude, perturbation 
theory is excellent in capturing the structure and evolu- 
tion of halos. By construction, the average value of the 
dynamical field is the same as the background field in 
perturbation theory. Hence, the modified drag term and 
modified mass terms are almost identical. Also, pertur- 
bation theory appears to be adequate for describing the 
evolution of the scalar field, so the effects of the dynamic 
field itself agree in the full theory. Additionally, since 
fluctuations dampen with time, we only need to satisfy 
perturbation theory constraints in the early universe. We 
see in Figure |4] that, for this particular amplitude of fluc- 
tuations, perturbation theory is excellent in describing 
the overall structure of the pancakes. 

Overall, as Figure|4]showed, perturbation theory seems 
to be adequate in describing the larger features of the 
pancakes: the number of caustics and the location of all 
three caustics at z = 0. With the notable exception of the 
fifth force, as we decrease initial scalar field amplitudes, 
perturbation theory becomes more accurate in providing 
an accurate solution. Even though the perturbation the- 
ory fifth force disagrees with the full theory by several 
percent, this is not enough to affect structure formation 
at our resolution. 

It is interesting to note that we can produce significant 
changes in structure formation even with a negligible fifth 
force, such as might happen with a sufficiently screened 
potential. Figure [5] shows the results from the full theory, 
but with no fifth force. In this case, we get stronger 
deviations in the innermost pancake substructure. Also, 
the second caustic forms slightly later than in standard 
cosmology, while the largest caustic is no different than 
when the fifth force is large. 

5. DISTINGUISHABILITY OF MODELS 

The models chosen in Table [T] are constructed to be 
consistent with current observations. For example, all 
these models predict values for the CMB peak locations 
that are within observational constraints. 

Figure [7] shows how each of our chosen potential pa- 
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Fig. 6. — Dark matter density profile at 2 = 0. The soUd line 
is the standard ACDM cosmology, and the dashed line is the final 
result in full theory when the fifth force is negligible. The numerical 
shot noise in the density near x = 1 is negligible and unimportant. 
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Fig. 7. — Evolution of (pi,, relative to the present-day value, for 
the various potential parameters labeled in Table [l] Note that the 
symbols are only to aid in distinguishing the lines. These calcula- 
tions are from perturbation theory. 

rameters in Table [T] affects the evolution of (p. These 
results agree with the perturbation results from Farrar 
and Peebles. There are in general two classes of viable (j) 
evolution tracks. We see from Table [T] that models such 
as A and C have a smaller initial value, but drop by 
roughly ten percent by the current epoch. On the other 
hand, models such as B and D have a high initial value 
and do not change much throughout cosmic history. 

These two classes of parameter sets have interesting 
consequences for the behaviors of the DM particle mass, 
the modified Hubble drag, the fifth force, and the dynam- 
ical field itself. Figure [8] shows how the mass ratio {ij) , 
the modified Hubble drag (7), the fifth-force acceleration 
iPpert), and the field energy density ratio {S) vary across 
the models in perturbation theory. For models of the sec- 
ond type, in which (f> does not change much, the mass at 
z = 50 is much closer to its present-day value, and hence 
this effect on structure formation will be smaller than 
in models where cf) varies strongly. Also, if (p does not 
vary much, then the Hubble drag will not be modified 
greatly. Together, this makes sense: a set of parameters 
that favors a slightly-varying scalar field will look more 
like standard ACDM than those which do not. Finally, 
for larger values of (p, the fifth force is much weaker, since 
the force is proportional to in perturbation theory. 

These behaviors play out accordingly in the final struc- 
ture at z = 0, as Figure [H] shows. These calculations were 
done using the full theory. When performing the analysis 
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Fig. 8. — Evolution of mass (rj), Hubble drag {7), fifth-force 
acceleration {Ppert)^ and dynamical field (<5) for the various poten- 
tial parameters labeled in Table [Tl The symbols are only to aid in 
distinguishing the lines. These calculations are from perturbation 
theory. 
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Fig. 9. — Dark matter density profile at z = 0. The solid line in 
all plots is the standard cosmology. The dashed line indicates the 
result in full theory, for the various potential parameters labeled 
in Table [1] 

in the full dynamical theory, we again have to be careful 
when setting the initial scalar field amplitude. For all 
cases, we set an initial contrast of 8^{z — 50) = 10~*. 

Models B and D do not differ much from standard 
cosmology or each other. Models A and C, while both 
different from ACDM, remain mostly indistinguishable, 
except that parameter set A, whose changes in </> are 
most drastic, produces caustics that are slightly farther 
away from the pancake midplane. Even though is in- 
creasing at the present epoch in parameter set C, this 
change happens at late times, and so it does not have a 
significant effect on the resulting structure. 

We also performed this study in perturbation theory, 
and we found that accounting for the effects of the full 
dynamics do not increase the distinguishability of the 
models. 

Surprisingly, we do notice some additional distinguish- 
ing features when examining the phase plots, as in Fig- 
ure [TOl Since models B and D remain nearly identical to 
concordance cosmology, we do not display them. While 
models A and C share many common features, the drag 
term 7 increases at late times in model C. This eventu- 
ally slows down particles at low redshifts and pushes the 
peak velocities closer to concordance cosmology. Hence, 
the higher peak particle velocities distinguish model A, 
whereas simply considering pancake density profiles may 
not. 
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Fig. 10. — Dark matter particle phase plots at z = 0. The solid 
line in all plots is from concordance cosmology, and the dotted lines 
are the results from various model parameters, which are labeled 
in the plots. These were calculated in the full theory. The symbols 
are only to aid in distinguishing the lines. 

6. CONCLUSIONS 

These one-dimensional simulations have clearly 
demonstrated that models of interacting dark matter 
and dark energy affect the growth and structure of 
plane parallel perturbations. Various consequences of 
these theories play important roles at different stages in 
halo evolution. Larger fifth forces and modified Hubble 
drag terms in the high-redshift universe greatly alter 
early structure formation, while an evolving particle 
mass can affect even late-forming structures. Ordinary 
quintessence, which is a dynamic field alone, does not 
significantly alter the formation of pancakes. Of all 
the effects, a modification to the Hubble drag gives the 
largest deviations in the resulting structure. 

We have found that models that are specifically tuned 
to match current observational constraints, such as the 
CMB peak locations and equation of state parameters, 
produce significant deviations from standard ACDM in 
the formation of pancakes. 

We have also found that an approach based on pertur- 
bation theory is adequate for understanding the general 
evolution of the scalar field. However, perturbation the- 
ory does not seem to be appropriate at high redshift. 
Even for arbitrarily small amplitudes, the fifth force in 
the full theory is larger than that obtained by perturba- 
tive techniques. This behavior deviates from the gen- 
eral pattern of perturbation methods, in that it does 
not accurately describe all characteristics for amplitudes 
smaller than a threshold value. 

This discrepancy arises from the fact that the fifth 
force directly depends on the gradient of the field. In 
the perturbative method described by Farrar and Pee- 
bles, the kinetic terms of the perturbation field equation 
of motion are dropped, producing a Poisson equation for 
the perturbation field and leading to the simplification 
of the fifth force expression found in equation (fT4| . This 
simplification is valid at low redshifts, but at high red- 
shifts the kinetic terms are still important. Thus, any 
realistic spectrum of perturbations for the scalar field 
which ignores the kinetic terms may be invalid for the 
very largest-scale perturbations. It also appears that the 
initial scalar field amplitude does not serve as a suitable 
"small" parameter for governing the appropriateness of 
perturbation theory. 

Because we can freely choose the initial amplitude, full 
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theory gives us more freedom to study the effects of vari- 
ous strengths of a fifth force. In the context of this model, 
this can also work in reverse: working within these frame- 
works, constraints on a fifth force in the dark sector could 
lead to limitations on the amplitudes of a DE scalar field 
in the early universe. 

Anninos and Norman have extensively studied 
the ef fects on baryons i n th e formation of pan - 
cakes (jAnninos fc NormanI 11994 lAnninos et al.l I1995D . 
While baryons appear to modify only the last-forming 
caustic, DM-DE interactions affect all caustics. Hence, 
the modifications to the last caustic may be masked by 
the baryon dynamics. Also, baryons do not appear to 
significantly modify dark matter particle velocities, while 
some of the models considered above produce significant 
variations. 

The structural differences described above are rela- 
tively small. However, in three dimensions these effects 
should appear as percent-level differences in the statisti- 
cal properties of large samples of dark matter halos. The 
various effects could produce different halo mass func- 
tions, which could be compared to other high-resolution 
m ass functions, such as the A CDM ones considered 
by I Warren et al.l (|2006( ). In fact, iMainini &: Bonomettol 
(|200l " have shown that DM-DE interactions do mod- 
ify analytic mass functions, and a direct simulation 
could be compared against these results. DM-DE in- 
teractions may also affect halo substruc ture, in which 
case they would impact NFW profiles (jNavarro et al.l 
Il997|) and dist ributions of conce ntrations in halo cata- 
logs (|Lukic et a l. 2007; Reed et a l. 2003). 

The present work focused on a single model and a sin- 
gle potential, simply to analyze the feasibility of a direct 
simulation approach. With a full three-dimensional sim- 
ulation, we can also examine the results of using other 
potentials, or even more complicated models. With fu- 
ture simulations, we can also contrast these structure 
formation results to the iV-body results from other the- 
ories of DE, such as m odified General R e lativi ty, which 
have been examined in lStabenau fc JainI (l2006l ). 

For the models with the largest deviations from 
ACDM, we notice that the sign of the difference of the 
caustic peak location changes. This implies that the 
shape and evolution of the mass function for these cos- 
mologies will change relative to ACDM, not simply the 
amplitude. This holds the most hope for observations, 
particularly for surveys such as the DES (jAnnis et al.l 



l2005f ). Most notably, we found that we will still get 
significant modifications to pancake structure even if 
the fifth force is negligible. Indeed, the modifica- 
tions to the mass function may even be more signif- 
icant without a fifth force. Since many recent ef- 
forts have concentrated mostly on cons t raining the fifth 
force (iKesden fc Kaniionkowskil l2006a| iFarrar fc RoserJ 
120071: iBertolami et al.ll2007t iGuo et aLlboOTl) . this result 
clearly shows that we should not ignore the other con- 
sequences, such as modifications to the Hubble drag, of 
these models. Also, the phase plots revealed another 
observable effect: the large-scale matter velocity field, 
which is also accessible to the DES survey. The feasibility 
of constraining these models with this behavior could be 
tested with velocity correlations in full three-dimensional 
simulations. 

However, there are several computational and theo- 
retical challenges in developing and analyzing full three 
dimensional simulations. For example, there should be 
a stronger theoretical understanding of the initial con- 
ditions of the field. While implicitly solving the scalar 
field is relatively straightforward in one dimension, re- 
quiring a simple Thomas algorithm, a full three dimen- 
sional solver would be much more complex, probably 
requiring a multigrid solver that could easily be com- 
bined with existing parallel a daptive- mesh co des such 
as FLASH (jFrvxell et al.l[2000t ) or GADGET-2 (|SDringe]l 
[200l . Since there might be interesting changes in halo 
substructure, simulations will require very high resolu- 
tion. Also, simulations will require many halos to get 
significant statistical results. 

Fortunately, none of these issues are intractable, and 
this approach holds much hope for providing a strong, 
consistent method of analyzing the nonlinear effects of 
these myriad dark energy models. 
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APPENDIX 

NUMERICALLY SOLVING THE SCALAR FIELD EQUATION 



We have 



dV p . 
d(f> (p 



a~^ = 0. 



(Al) 



We will let V"(0) — Kcj) and we have re- introduced c for clarity. We may break up equation I|A1|) using standard 
operator splitting techniques: 







(A2) 
(A3) 



Here and below, superscripts are temporal indices and subscripts will denote the index of the location on the mesh. 



We solve Eq. (|A2[) by reducing it to a set of two Ist-order equations: 



(A4) 



dt 

%r—-^'r- (A5) 

Since the sound speed is so high (vs = c/a), we found that a direct integration scheme was highly unstable except 
for unworkably small At. So, we use a midpoint method for Eq. (|A4[) and a fully implicit method for equation (|A5|) : 



i = 0r + ^At(0r + 0r^) (A6) 
^(0i^\-2</>«+<^W). (A7) 



We substitute our expression for (f>l into equation (|A7|) . yielding the matrix equation Ax — 6, where 

X^ = ^P, (A8) 



and 



A 



6, = 0r + a0r+i+7f^'/'r+i 



1 + 2(77 — cr7 • • • — CT7 \ 
—(77 1 + 2(77 — <77 ■ • • 
—(77 1 + 2(77— (77 



(A9) 



\ -(77 

Note the periodic boundary conditions. We have defined 



-(771 + 2(77/ 



At 



(AlO) 



We solved this matrix- vector equation exactly u sing a Thornas alg orithm (jConte k deBooi]|1972[ ) , modified for periodic 
boundaries via the Sherman-Morrison formula (jPress et al.lll993) . 
We may now update our solution through equation (jA3p by expanding in a Taylor series: 



WAi + i^At2 + 0(Ai3), 



where 



and 



dt^ 



-12 



+ 3-^(a")-(a"). 



(All) 



Here, p = (p" - /dt. 
At any step ri, after computing (/>" across the grid, we may update the scale factor using a Taylor series: 



a" + d"At + -a"At2 + 0{At^), 



(A12) 
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where a and a are given by the Friedmann equation and Friedmann energy equation, respectively: 

1/2 



a = 



SttG 



"0 



-n 1 2r. -2 47rG' 

a =--HQilm^a 5— 

Ignoring fluctuations, the equation of state parameter, w, is 

^ -r2 

P4> 



(1 + 2,w)p^a. 



W 



P4- 



(A13) 
(AM) 

(A15) 



REFERENCES 



Alcaniz, J. 2006, Braz. J.Phys., 36, 1109, astro-ph/0608631 
Amcndola, L. 2000, Phys. Rev. D, 62, 043511 

Amcndola, L., Campos, G. C, & Roscnfcld, R. 2007, Phys.Rcv.D, 
75, 083506 

Amcndola, L. & Tocchini-Valcntini, D. 2001, Phys. Rev. D, 64, 
043509 

Anninos, W. Y. & Norman, M. J. 1994, ApJ, 429, 434 

Anninos, W. Y., Norman, M. L., & Anninos, P. 1995, ApJ, 450, 1 

Annis, J. et al. 2005, astro-ph/0510195 

Bean, R., Carroll, S., & Trodden, M. 2005, astro-ph/0510059 
Bertolami, O., Pedro, F. G., & Delliou, M. L. 2007, astro- 
ph/0705.3118 

Conte, S. & deBoor, C. 1972, Elementary Numerical Analysis 

(McGraw-Hill, New York) 
Corasaniti, P.-S., Hutercr, D., & Melchiorri, A. 2007, ArXiv 

Astrophysics e-prints 
de Bernardis, P. et al. 2000, Nature, 404, 955 
Dodelson, S., Kaplinghat, M., & Stewart, E. 2000, Phys.Rev.Lett., 

85, 5276 

Doran, M. & Jackel, J. 2002, Phys. Rev. D, 66, 043519 
Farrar, G. R. & Peebles, P. 2004, Astrophys.J., 604, 1, astro- 
ph/0307213 

Farrar, G. R. & Rosen, R. A. 2007, Phys.Rev.Lett., 98, 171302 

Fryxell, B. et al. 2000, Astrophys.J. Supp, 131, 273 

Gnedin, N. 1993, in Bulletin of the American Astronomical Society, 

Vol. 25, Bulletin of the American Astronomical Society, 1464 — h 
Guo, Z.-K., Ohta, N., & Tsujikawa, S. 2007, Phys. Rev. D, 76, 

023508 

Hockney, R. W. & Eastwood, J. W. 1988, Computer simulation 

using particles (Bristol: Hilgcr) 
Hucy, G. & Wandclt, B. D. 2006, Phys.Rcv.D, 74, 023519 
Ivezic, Z. 2007, in American Astronomical Society Meeting 

Abstracts, Vol. 210, American Astronomical Society Meeting 

Abstracts, 66.05 

Kesdcn, M. & Kamionkowski, M. 2006a, Phys.Rev.Lett., 97, 

131303, astro- ph/0606566 
— . 2006b, Phys.Rcv.D, 74, 083007, astro-ph/0608095 
Koivisto, T. 2005, Phys.Rcv.D, 72, 043516 

Lukic, Z., Hcitniann, K., Habib, S., Bashinsky, S., & Ricker, P. M. 
2007, ArXiv Astrophysics e-prints 



Maccio, A. V., QucrccUini, C, Mainini, R., Amcndola, L., & 
Bonometto, S. A. 2004, Phys. Rev. D, 69, 123516 

Mainini, R. & Bonometto, S. 2006. Phys. rev. D, 74, 043504, astro- 
ph/0605621 

Mancra, M. & Mota, D. F. 2006, MNHAS, 371, 1373 
More, J. G. 1988, PhD thesis, AA(Edinburgh Univ. (Scotland).) 
Navarro, J., Frcnk, C, & White, S. 1997, Astrophys.J., 490, 493 
Nusscr, A., Gubscr, S., & Peebles, P. 2005, Phys.Rev.D, 71, 083505, 

astro-ph/0412586 
Olivares, G., Atrio-Barandela, F., & Pavon, D. 2006, Phys.Rev.D, 

74, 043521 

Peebles, P. & Ratra, B. 1988, Astrophys. J.Lett., 325, L17 
Percival, W. J. et al. 2001, MNRAS, 327, 1297 
Perlmutter, S. et al. 1999, ApJ, 517, 565 

Press, W., Flannery, D., Teukolsky, S., & Vcttcrling, W. 1992, 
Numerical Recipes in C: The Art of Scientific Computing 
(Cambridge University Press) 

Ratra, B. & Peebles, P. 1988, Phys.Rev.D, 37, 3406 

Reed, D., Gardner, J., Quinn, T., Stadel, J., Fardal, M., Lake, G., 
& Governato, F. 2003, MNRAS, 346, 565 

Riess, A. G. et al. 1998, AJ, 116, 1009 

Springel, V. 2005, MNRAS, 364, 1105 

Stabenau, H. & Jain, B. 2006, Phys.Rev.D, 74, 084007, astro- 
ph/0604038 

Valinia, A., Shapiro, P. R., Martel, H., Vishniac, E. T., & 
Villumsen, J. V. 1995, in Bulletin of the American Astronomical 
Society, Vol. 27, Bulletin of the American Astronomical Society, 
1365 

Warren, M., Abazajian, K., Holz, D., & Teodoro, L. 2006, 

Astrophys.J., 646, 881, astro-ph/0506395 
Yuan, W., CentrcUa, J. M., & Norman, M. L. 1991, 

Astrophys.J. Lett., 376, L29 
Zel'dovich, Y. B. 1970, A&A, 5, 84 

Zimdahl, W. & Pavon, D. 2001, Phys. Lett. B, 521, 133, astro- 
ph/0105479 



